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It has been shown recently that the amplitude of the creation of n real scalar 
particles by one virtual boson near n-particle threshold exhibits exponential 



behavior at n ~ 1/A. We extend this result to the processes of multiparti- 



cle production at threshold by two virtual bosons. We find that both the 
tree-level amplitude and leading-ra loop corrections have the same exponen- 
tial behavior, with the common exponent for 1 — > n and 2 — > n processes, 
for various kinematics. This result strongly indicates that the exponent for 
multiparticle amplitudes is independent of the initial state and that there may 
exist a semiclassical approach to the study of multiparticle production. 
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I. INTRODUCTION 



A reliable calculation of multiparticle amplitudes remains an interesting problem in quan- 
tum field theory. The problem exists in most weakly coupled theories including the 
simplest case of ^<£> 4 model, where it has been observed that the tree amplitude to produce 
n final particles from an initial virtual one exhibits a factorial dependence on the number 
of outgoing particles like n!A n / 2 , which leads to an unacceptably rapid growth of the tran- 
sition rate at the tree level at large enough n, n ~ 1/A. Since the contribution from the 
first loop is of order Xn 2 as compared to the tree-level result, and from the k-th loop one 
expects a contribution of order (Xn 2 ) k , it is clear that the calculation of the amplitude at 
n ~ 1/A requires the summation of the whole perturbation series. Recently, it has been 
demonstrated that at the n-particle threshold, leading in n contributions from each loop 
sum up to exponent [[K]] 

A^ n = A^ n -e BXn2 (1) 
where B is some numerical constant and 

ATJn = n\ (~J 2 (2) 

is the tree amplitude (the scalar boson mass is set equal to 1). Equation (P determines 
correctly the threshold amplitude with the account of all loops at n ~ 1/ V% in which case 
the loops have the same order of magnitude as the tree-level contribution. More importantly, 
at n ~ 1/A, Eq. ([[]) provides strong indication for the exponentiation of the loop corrections. 
In other words, one expects that at n ~ 1/A, the threshold amplitude has the following form, 

A(l -> n) v^! • e^ F(An) (3) 

where F is some function that has the following expansion at small An, 

F(Xn) = F tree + BX 2 n 2 + 0(A 3 n 3 ) 

with 
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so at An <S 1, Eq. @ reduces to Eq. (p. The behavior of the function F is presently 
unknown at An ~ 1. 

It is likely that the exponential form of the amplitude, Eq. @, is a consequence of 
the semiclassical nature of the processes of multiparticle production. So, one may hope 
that there exists some semiclassical-type technique for calculating the amplitudes of these 
processes. 

In this direction, one interesting suggestion 0,|11|,0 is to try to generalize the Landau 
method for the calculation of semiclassical matrix elements from one-dimensional quantum 



mechanics to field theory. The Landau method |T3[ is a powerful technique which allows one 
to calculate the matrix elements of almost any regular operator between two semiclassical 
states, with different energies, of a particle moving in one-dimensional potential. In field 
theory, the direct analogues of these matrix elements are the multiparticle amplitudes: for 
instance, the amplitude to produce n scalar bosons at threshold by a virtual particle can be 
written in the form (n|0|O), where |n) is the state with n particles at threshold. 

One specific feature of the semiclassical matrix elements in one-dimensional quantum 
mechanics which is captured by the Landau method is that they do not depend, to the 
exponential precision, on the operator for which they are calculated; they depend only 
on the states between which the operator is sandwiched, as well as on the details of the 
potential. Therefore, if there really exists a generalization of the Landau method to field 
theory, one would expect that the exponents of the matrix elements (n|O|0) should not 
depend on the operator O, provided the latter does not depend on the coupling constant 
A. In other words, the multiparticle amplitudes should be the same, to the exponential 
precision, for all few-particle initial states. So, a check of the existence of the Landau-type 
procedure for calculating multiparticle amplitudes will be the demonstration that they are 
indeed independent of the initial state. Note that similar conjecture has been made for 



instanton-induced processes [14,15 



Presently, it is unclear how to perform reliable calculations in the regime Xn ~ 1. How- 
ever, the technique of Ref. |I0] can be generalized to deal with initial states that may contain 
more than one virtual particle. In this way one should be able to verify that the exponential 
part of the tree level amplitudes is the same as in Eq. (|]) for all few-particle (containing much 
less than 1/A particles) initial states, and that the leading-n loop corrections exponentiate 
to exp(BXn 2 ) with exactly the same coefficient B as in Eq. (p]). 

In this paper we perform this check for initial states that contain two virtual bosons and 
final states containing n real bosons at rest (p/ = 0). Since there are two incoming particles, 
the energy can be distributed arbitrarily between them. The first case that we consider is 
the amplitude 2 — > n integrated over 4-momentum of one initial particle. This quantity is 
equal to the matrix element (n\ip 2 \0) at n-particle threshold. In two other cases we evaluate 
the amplitude of scattering of two particles with 4-momenta (E, p) and (n — E, — p) into 
n particles at rest (we set the mass of the boson equal to 1) in two different regimes: at 
E <C n, so that the energy of one initial particle is much smaller than that of the other, and 
at E of order n when the energy of the two initial particles are comparable to each other 
(we will take the spatial momentum p of the incoming particle to be small, |p| ~ 1). In 
the latter case we consider the regime when E does not belong to the interval (0, n), so one 
initial particle carries negative energy (it may appear more natural to view this particle as 
an outgoing one; this is of course a matter of terminology: this particle is virtual anyway, 
and we are free to call it incoming). The interval < E < n is peculiar: for these values 
of E even the tree-level amplitude has singular behavior as a function of E, which is a 
consequence of the threshold kinematics, and does not have a well defined limit in the large 
n, fixed E/n regime. Due to this peculiarity, the case < E < n is not suitable for our 
main purpose. 

The technique that we will apply is a direct extension of one developed in Ref. |TIJ. In 
all three cases, the calculations show that the tree level amplitude and the leading-ra loop 
corrections exponentiate to the same exponent as in the case of 1 — > n process. This result 
strongly indicates that the exponent of multiparticle amplitudes in fact do not depend on 

4 



the initial state and supports the existence of some semiclassical, possibly Landau-type, 
approach to the study of multiparticle production. 

The paper is organized as follows. In Section [TI| we will briefly review the approach 
developed in Ref. ]ID| for summing leading-ra loop corrections. In Section [ITT], this technique 
is applied for calculating matrix elements of the operators tp 2 . Section [TV] is devoted to 
2 — > n amplitudes in the case of one soft initial particle, while in Section [V| the case of two 
hard initial particles is considered. Section [Vl| contains concluding remarks. Appendix is 
devoted to a brief discussion of the case < E < n, and the tree and one-loop results are 
presented. 



II. GENERAL FORMALISM 

We consider <£> 4 theory with unbroken discrete symmetry in (d+1) dimensional spacetime 
with the Lagrangian 

L = \{d, v f - - ^ 4 (5) 

Hereafter we set the mass of the boson m to be equal to 1. 

For calculating multiparticle amplitudes at threshold, there exists a convenient formalism 
0] which reduces the problem to the calculation of Feynman graphs in certain classical 
background. In Ref. [|j] this formalism has been developed for the 1 — > n amplitude, but it 
is easy to extend it for treating processes with two incoming particles. 

Let us outline briefly this technique (see Ref. |4| for details). Consider a transition from 
an initial virtual particle with (d + l)-momentum = (n, 0) into n final particles, each 
with (d + l)-momentum (1,0). The reduction formula for the amplitude can be written in 
the following form, 

•Xfl 

(n|0(x)|O> = lim lim (eg - 1)"-— (0|^(x)|0> J=Joexp(ieot) (6) 

£q^1 Jo— »0 OJq 

where the matrix element is calculated in the presence of a source J = Joexp (ieot). 
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Taking the limits Eq — > 0, Jo — > simultaneously, one can show that A n is determined 

by 

An = Q^mo)\ z=0 (7) 

where the expectation value (0|</?|0) is calculated in the following classical background 

which is a solution to the field equation. 

It is straightforward to generalize this technique to the case of matrix elements of arbi- 
trary operators. In particular, the matrix element of the operator (p 2 can be obtained by 
differentiating its vacuum expectation value in the presence of the background field varphi , 

H^|0> = ^<% 2 |°> 1,0=0 (9) 



Analogously, to calculate the amplitude of scattering of two initial particles with momenta 
(E, p) and (n — E, — p) into n bosons at threshold, one should differentiate the corresponding 
full propagator n times with respect to z , 

A 2 ^ n {E, p) = -^-j d d+l x d d+l y e * E « +*(«-B)v e<p(y-*)x)(ar, y) (10) 

where T>(x, y) is the two-point Green function calculated in the classical background (pQ. It 
is convenient to use the mixed coordinate-momentum representation, 



V p {x°,y°) = Jd d xd d ye l ^-^V(x,y) 



so Eq. (|T0"D can be rewritten in the following form, 



A 2 ^ n (E,p) = Jdx dy e* E * 0+ ^y°V p (x .y ). (11) 

So, the amplitudes to produce n final particles at threshold can be obtained by differenti- 
ating the corresponding Green functions ("generating functions"), calculated in the presence 
of the background field <po, n times with respect to the parameter of the background, zq. 
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The Green functions which enter the right hand sides of Eqs. (j^j and (|TT|), in their turn, can 
be computed in the perturbation theory around the classical background <p . It is convenient 
to introduce the Euclidean time variable, 

1, A , 17T 

T = it + -\n- + \nz + —, (12) 
in terms of which the background field has the form 

Mr) = - ? /!odh? 

Note that the background field has a singularity at r = in/2. By expanding the Lagrangian 
around the background (po one obtains the Feynman rules shown in Fig. |I[ 

Obviously, perturbative calculations become more and more complicated at higher loops. 
However, if one is interested only in the leading-n contribution from each loop level, consid- 
erable simplification occurs ||10|| , which allows for the summation of the whole perturbative 
series. The key point is that, at each loop level, the large-n behavior of the amplitude 
depends only on the structure of the singularity at r = in/ 2 of the generating functions (in 
our case (0|v? 2 |0) or V p (x°,y )). Thus, our strategy is to find the generating function near 
the singularity and after that recover the multiparticle amplitude. 

In the case of 1 — > n process, the procedure has been developed for obtaining the 
leading singularity of the corresponding generating function at r = in/ 2 at any given loop 
order, which reduces the problem to the calculation of tree graphs in some effective theory. 
We will apply this technique to the case of the processes with two initial particles. Let us 
begin with the matrix element of (p 2 . 



III. MATRIX ELEMENTS OF (p 2 

Let us first discuss the Feynman graphs that contribute to the generating function, 
(0|<£> 2 |0). There is one tree-level graph shown in Fig. 0. At the one-loop level, there are two 
different types of graphs. The two-loop graphs are much more numerous, only one of them 
is presented in Fig. [|]a. 
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As explained above, at each loop level it is sufficient to calculate only the leading sin- 
gularity of the generating function. The latter can be found by extending the technique of 
Ref. ||10|, that reduces the problem to the evaluation of tree graphs. Since this extension is 



straightforward, let us summarize here only the prescription (see Ref. |10| for details): 



• For a given loop graph one should cut some propagator lines in such a way that the 
resulting graph is tree and connected. In cases when there are various ways to cut, 
the result is given by the sum over all possibilities. The number of lines to be cut is 
equal to the number of loops. 

• For each propagator line, say, G(t, t'), that has been cut, one attributes the factors 
B l l 2 \Lp 2 {r) and B 1 / 2 \(p 2 (r') to the two vertices that this line connects, i.e., at r and 
t', respectively. The constant B is defined by 

o = f^JP 9 

J (2tt)« 8u p (ul-l)(ul-4) 

where u p = \/p 2 + 1. All propagators that have not been cut should be replaced 
by the operator (d 2 + 3\ipl)~ 1 . Therefore, the tree graphs obtained by our cutting 
procedure do not contain spatial momenta explicitly. 

Let us demonstrate how this prescription works for a particular graph shown in Fig. 
|3]a. This graph has a symmetry factor of 1/4. Since it is a two-loop graph, the number of 
propagators to be cut is 2. There are in fact 5 possibilities to cut in such a way that the 
graph remains connected: one can cut the lines (1,3), (2,3), (1,4), (2,4), or (3,4) (one cannot 
cut two propagators (1,2), since then the graph will become disconnected). In fact, one can 
see that the first four ways lead to essentially the same graphs, while the graph obtained 
in the fifth case has another topology. As a result, one obtains a sum of two graphs with 
symmetry factors of 1 and 1/4, as shown in Fig. ^b. The black circles ("bullets") represent 
the factors B x l 2 \kp 2 . 

By analyzing the tree graphs obtained by cutting the initial loop graphs, Figs. |2] and |3|, 
one can see that they have the form of a tree cascade starting from two initial lines corre- 
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sponding to the two operators <p(r) ( a two-branch tree). Each branch contains contributions 
of <£>o, B l l 2 \ip^ and higher power of A. One observes (analogously to the 1 — > n case) that 
the same series of graphs is obtained when one solves the field equation without the mass 
term, 

<9^ + Ay? 3 = (13) 

perturbatively with respect to A, with the condition that the first two terms in the expansion 
over A are (fo + B l l 2 \(pl. 

Eq. (p~3| ) has the following exact solution 

^ = -jf— (14) 

V A T - T 

which has a single pole at r = tq. One can fix To by requiring that the solution has the 
expansion ip c \ = y? + B 1 / 2 \ip' 2 l + ■ ■ ■. As a result, one obtains 



r = wr/2 - V2BX. 
Then ip c \ can be represented in the following form, 

^ cl = ^oE( 5l/2 Vo) fc (15) 

k 

At first sight, the generating function is equal to (p 2 ^ which has the following expansion, 

oo 

^ = ^£(* + i)(a 1/a W' (is) 

fc=0 

However, this is not true. First, one notes that all graphs obtained by the cutting procedure 
contain even number of bullets, so there should be only even powers of B l l 2 X in the series for 
the generating function. Second, from every tree graph with 21 bullets one can reconstruct 
(21) ! / (2'/!) graphs of the original theory at /-loop order by pairing the bullets into propagator 
lines. So, to recover the generating function one should omit in Eq. (|IED all terms with odd 
k and for terms with even one should multiply the coefficient by the factor 

(2/)!/(2'/!). In this way one obtains, 



<(V|o> = ^£ 



1=0 



oo 



(2Z+1)! 
2 Z Z! 



(X 2 B^) k 



(17) 



To obtain the matrix element {n\ip 2 \0) one substitutes Eq. (pi) into Eq. ([H]) and differentiates 
n times with respect to z. Recalling Eq. fll2"D , one obtains the following result, 



can be identified with the matrix element calculated at the tree level, which differs from the 
tree 1 — > n amplitude @ only by a pre-exponential factor proportional to n. The leading-n 
loop contributions exponentiate to the same factor of exp(BXn 2 ). Therefore, the exponent 
for the matrix element of (p 2 coincides with the 1 — ► n amplitude, which is the desired result. 

IV. AMPLITUDES 2 -> n WITH ONE SOFT INITIAL PARTICLE 

In this section we consider the process of scattering of two initial virtual particles, among 
which one is soft (with energy E <n) and the other is hard, into n final particles at threshold 
(the spatial momenta of initial particles are assumed to be small, i.e. of order 1). 

It is convenient to rewrite Eq. (]TT|) in a slightly modified form. One notes that it is suffi- 
cient to consider only one integral over dy°, the integral over dx results in the delta-function 
of energy conservation, which we will not explicitly write in what follows. Furthermore, one 
introduces the notation z{t) = zoe lt , and goes to Euclidean time, then the 2 — > n amplitude 
is written as follows, 




where 




G(t) 



dz 



where 




(18) 
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We will call G(r) the generating function for the 2 — > n amplitude. At the tree level 
this result was obtained firstly in ]TB| by direct summation of graphs. Let us begin with 
discussing the tree level. 



A. Tree level 

At the tree level, it is easier to find the amplitude and then recover the generating 
function. The tree-level amplitude can be found directly from the exact tree propagator in 
the background field (p , Eq. (g) [0,0 

D p (t,t') = ^rWWftWeV - r) + f^r)fnr')e(r - r')) 

where 



W p = 2w(u 2 - l)(u 2 -4), u = Jp 2 + l 



and 



(t) = e^ 



12 6(w-2) , 

^V2 + TTe^ + {uJ ' l){uJ ' 2) 



= /fir) 



To evaluate the integral (fTH), one expands the functions fi and f 2 in series, 



/r(r) = e"^(-l) fc /i fc e 2fa 

fc=0 



where 



/ 2 w (r) = e- WT V(-l) fc / 2fc e 2fcr 



fc=0 



/i* = - l)(w - 2) + 6cu + 12fc 



/ 2fe = + l)(w + 2) - 6cu + 12A; 



(19) 
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The tree amplitude 2 — > n can be now obtained by direct calculation. One finds, 

Atree ( F , _ , ( X Y /2 1 flkf^-h) hkfl(^-k) \ 

It is worth noting that A^ n (E, p) develops a series of poles at E = 2k ± to. The physical 
reason for these poles is that for these values of E the tree graphs may include a propagator 
with on-shell momentum (see Fig. |4j). 

To obtain the asymptotic behavior of the amplitude at large n, one notes that at E ~ 1, 
the sum in Eq. (|20|) is saturated by a finite number of terms with k ~ 1. So we can replace 
/i,2(n/2-Jfe) by 6n and extend the sum to = oo. We obtain, 

/ \ \ n /2 i oo / f f \ 

Itree f tp \ c„l„ I A 1 1 V^/ Jl* J2k \ 



A«Z(E,p) = 6n\nl-\ — £ 



Wp^V 2fc + c,; - E 2k-uo-E) 
Note that the sum in this equation converges and does not depend on n. Let us introduce 
the function, 

no? „\ 3 1 v' ^ lfc ^ 2fc i 



2 W p \2k + u - E 2k-u-E i 



The tree amplitude can now be rewritten in the form 

A^ n (E,p) = 4n-n\(±) ' C(E,p) (21) 



Note that the tree amplitude is proportional to the tree matrix element of ip 2 , 

Af™ n (E,p) = \C(E,p)(n\p 2 \0) tree 

and correspondingly differs from the tree 1 — > n amplitude (|2|) only by a factor of n. There- 
fore the generating function is proportional to that of the operator ip 2 , 

G trce (r) = AC(£,p)^( r ) (22) 

One can verify by differentiation that the generating function (^) indeed gives rise to the 
amplitudes (|2l|) at large n. 
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B. Loop corrections 



To sum leading-n loop corrections to the 2 — > n amplitude, one should calculate the 
generating function 

G(t) = e Er Jdr'V(T,T')e- ET ' 

in the region close to singularity. The technique that we apply here is the same as that 
used in the case of the operator (p 2 : one cuts every loop graph so that it becomes tree 
and connected and attaches bullets to the cut lines. In this way one obtains the leading 
singularity of the propagator T>(t, t') as a sum of graphs that typically have the form shown 
in Fig. H One can see that these graphs are the same as those obtained by expanding the 
propagator in the background field <p c \ 

D c1 (t,t') = (-dl + 1 + 3A0C!)" 1 

in the perturbation series in A, where the background field <p c \ and its expansion are given by 
Eqs. ([L4]) and (|T5|), respectively. The only difference is that only graphs with even number 
of bullets is present for T> and the symmetry coefficients differ by a factor of (2l)\/(2 l l\), 
where 21 is the number of bullets. 
The quantity 

G c1 (t) = e ET J dr'D cl (T, r')^' (23) 

can be easily calculated since the fields (p c \ and (fo have the same behavior around their 
singularities, while the singularities are located at different points. Recalling Eq. (|22"D one 
writes 

G c i(r) = \C(E, P )<&(t) 

oo 

= \C(E,p)rtJ2(k+l)(\ 2 B V l) k (24) 

1=0 
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G(t) can be obtained from this expression by dropping all terms in the sum with odd k and 
correcting the coefficients by the factors (2Z)!/(2'Z!). Comparing with Eq. (|ITD, one finds 
that the generating function in our case is proportional to that of the matrix element of ip 2 , 

G(t) = AC(£,p)((V|0> 

so the loop corrections exponentiate in the same way as for (n\(p 2 \0), 

A 2 ^ n (E,p) = A^ e n (E,p)e BXn2 

So, we have established that both the tree expression and the leading-n loop corrections 
for the amplitude 2 — > n in the case when one initial particle is soft sum up to exp(^F tree + 
BXn 2 ). 

V. AMPLITUDES 2 -> n WITH TWO HARD INITIAL PARTICLES 

In this section we consider the case when both initial particles have energies of order n. 
For the amplitude to have a regular limit at large n, we take the energy of one incoming 
particle, E, to be negative, while the energy of the other is larger than n. The case when 
energies of both initial particles are positive is briefly discussed in Appendix. 

A. Tree level 

The tree amplitude in our case can be determined from its representation in terms of the 
sum in Eq. fl2"0p. However, from Eq. (^0) it is not straightforward to extract its asymptotic 
behavior in the regime n — > oo, E/n = fixed. We adopt here another approach. We will 
find the generating function from the equation that it obeys, 

(-d 2 + 3Atfj) (e- ET G trcc {r)) = e~ Er (25) 

(we have omitted the term uj 2 which is inessential near the singularity). Since the amplitude 
to produce n particles depends on the details of the behavior of generating function in the 
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region |t — m/2\ ~ 1/n, and E ~ n, we look for the solution of Eq. (|25|) in the region 
\r — in/2\ ~ 1/.E. Near the singularity, Eq. fl25|) reduces to 



H + (7db^)( e " TG "" (r ^ e " T (26) 

Eq. ( PB| ) can be solved exactly, the solution has the following form 

G "~M = "BP ( P(T - ii;2)i + 1(7^) + 3 + B(T - " /2) - B2(T - n ' 2f - 

-E^r-i.^'j^j-A (27) 

— OO / 

We have kept in Eq. (^) not only the leading singular term, but all terms which are of 
order E~ 2 at E\t — 7r/2| ~ 1. Those terms are equally important for the evaluation of the 
amplitude: while the leading singular term (the first term in parenthesis in Eq. ([27])) gives 
rise to the contribution of order nE~ A , the second term produces E~ 3 , etc. (see below). 
Note that E < 0, so the integral in this equation converges. 

The large-n asymptotics of the amplitude is uniquely determined by the behavior of the 
generating function which has been found in Eq. (|27l). To recover the amplitude one should 
find a series in e 2r which has the same behavior around the singularity as the r.h.s. of Eq. 
(p7|). Apparently there are various ways to write down the series in e 2r that coincides with 
Eq. ([Z7]) near the singularity, but they all yield the same result for the amplitude. Technically, 
the most convenient way to write the series is to make the following replacement, 

1 + e 2r 

T - m/2 -> _ (28) 

which is valid near the singularity, so the generating function obtains the form 



E 3 r p~ Et ' 
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which is indeed a series in e 2r . Recalling Eq. (0), the tree amplitude can be derived by 
differentiating with respect to z. One obtains, 

A ^- l ^fm^W^-E? + E>) ,29) 

Again, the exponential part of the amplitude is equal to exp(A _1 F tree (An)), while the energy 
dependence enters only the pre-exponential factor. 



B. Loop corrections 

In analogy to the case of one soft initial particle, one should first calculate G c \, Eq. (^3|). 
To do this we recall that D c \ is the propagator in the background field f c \ which has the 



pole at r = in/2 — y2\B. So, instead of Eq. ( p6|) one has 



-dl + - m „ ) (e-^GUr) = e~ Er 



6 

(r-m/2 + V2XB) 2 , 

Since this equation has the same form as Eq. (|26|), its solution is given by a formula analogous 



to Eq. (P7|), with the only difference that the pole at iir/2 is shifted to m/2 — \J2XB. 
Making the replacement fl28|) , one obtains the solution in the following form (we do not 
write explicitly the terms that are regular at the singularity) 
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5E2 U 2 fl + e 2 --2v / 2TB) E (l + e 2 - - 2V2XB 



G(r) 



E 3 3 } e~ El ~' \ 

(l + e 2r - 2V2XB) e Er ' / t dr' + regular terms (30) 



which should be differentiated with respect to z in order to obtain the amplitude. Remem- 
bering that we should omit terms with odd powers of B x l 2 and multiply the coefficient of 
(B 1 / 2 ) 21 by (2/)/(2'Z!), we obtain after a tedious but straightforward calculation, 

1 tree (XBn 2 ) k _ tree \Bn 2 



A _ /(tree y/v-^'t. > _ A 



fc=0 



2-m 



C 
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where A^ n is given by Eq. fl29|). Therefore, although in the case of two hard initial particles 
the calculations are more complicated, the result is the same as in the case of the matrix 
element of if 2 or when one initial particle is soft: leading-n corrections sum up to the 
exponent exp(BXn 2 ). 

VI. CONCLUSION 

We have studied the amplitude of the processes 2 — > n in the <£> 4 theory by the technique 
that allows to sum up all leading-n loop corrections at n-particle threshold. We have shown 
that the 2 —>■ n amplitudes, regardless of how the initial energy is distributed among the 
two initial particles (except for some peculiar cases), coincide with the amplitude 1 — > n 
to the exponential precision, at least when only leading-n loop contributions are taken into 
account. The similar result can be easily obtained in the case of broken discrete symmetry. 
Our results, though not being a rigorous proof, strongly support the hypothesis that the 
amplitude of multiparticle processes is independent, in the exponential approximation, of the 
initial few-particle state. The picture here is similar to that of quantum mechanics, where 
the calculation of semiclassical matrix elements by the Landau method requires no knowledge 
on the details of the operator sandwiched between the semiclassical states. Our calculations, 
therefore, indicate that there may exist an extension of the Landau method to field theory, 
which, hopefully, could bring about the understanding of multiparticle amplitudes. 

ACKNOWLEDGMENTS 

The authors are indebted V. A. Rubakov for numerous helpful discussions and stimulating 
activity. We thank Yu. M. Makeenko and P. G. Tinyakov for valuable discussions. D.T.S. 
thanks Rutgers University for hospitality. This work is supported in part by INTAS grant 
# INTAS-93-1630. The work of D.T.S. is supported in part by Russian Foundation for 
Fundamental Research, grant # 93-02-3812. 



17 



APPENDIX A: 2 n AMPLITUDES FOR < £ < n 



In this Appendix we will briefly consider the behavior of the tree amplitude 2 — > n in 
the case when both initial particles have positive energies of order n, and also present the 
result for the one-loop correction to this amplitude. 

The tree amplitude can be derived by making use of Eq. (|2"0|). At E ~ n and < E < n, 
the sum in Eq. ( PU| ) is saturated by the terms with k ~ E/2. For these values of k one has 
fi,2k ~ 6E, /i,2(n/2-fe) ~ Q(n — E). The sum in Eq. ( pOl) can be extended so that k runs from 
— oo to oo, and one writes, 

' X\ n/2 36E{n - E) ~ / 1 1 



k^yZk + u - E 2k-u-E 
, /A\ n/2 72E(n-E) ( it , ^ , vr /rn \ 

Note that the amplitude has a very rapid oscillating behavior as a function of E, so it has 
no regular limit in the regime n — > oo, J5/n = fixed, but rather the amplitude depends on, 
say, the fractional part of £7/2. Due to this behavior of the amplitude at the tree level, one 
should not, in general, expect that the loop corrections will sum up into a regular factor like 

e BXn 2 

We have performed the calculation of the one-loop correction to the amplitude by direct 
evaluation of the three graphs shown in Fig. H. At large n, the result is 



/ \ \ ™/ 2 io 

A i-ioo P = n! I ^ J —XE(n - E) 



[e 2 + (n- Ef) (cot \{E + u)- cot |(£ - uj) + 



+ 9E(n -E) J (cot |(£ + cu k + cu p _ k ) - cot - cu k + w p _ k )) 



(Al) 



One sees that the one-loop correction is rather complicated and in general is not equal to 
BXn 2 ■ A tree , unlike the cases considered in the body of this paper. The only exception is 
the theory in (2+1) dimensions, d = 2. In that case the integral in Eq. (|A1[ ) is infrared 
divergent at k = and k = p, the same is true for B. If one introduces an infrared cutoff po 



(say, by considering final particles not exactly at threshold but with finite momenta of order 
Po), Eq. reduces to y 4 1_loop = BXn 2 ■ A tree . This result can be easily understood since in 
(2+1) dimensions, the loop corrections near threshold are dominated by the rescattering of 
final particles. The analysis then is completely similar to the case of 1 — > n processes, and 
further details can be found in Ref. JTD|. 
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FIG. 1. Feynman rules for calculating multiparticle amplitudes 





FIG. 2. The tree and one-loop graphs which contribute to the matrix element (0|<^ 2 |0). 
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FIG. 3. The two-loop graph for the operator ip 2 (a) and a representation (b) of its leading 
singularity 
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FIG. 4. A typical tree graph for the process 2 — ► n (ordinary Feynmann rules are assumed). 
When the momentum running along one of the thick lines is on-shell, the tree amplitude is infinite. 




FIG. 5. A typical graph in the perturbative expansion of P(r, r') 




FIG. 6. The one-loop graphs contributing to V(x,y) 
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